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DETERMINATION  OF  PERFORMANCE  PARAMETERS  FOR 
FIN-STABILIZED  FREE-FLIGHT  MISSILES 


ABSTRACT 


The  methods  used  by  the  Firing  Tables  Branch  for  the  analysis 
of  flight  tests  of  fin-stabilized,  free-flight  missiles  are  outlined. 
Examples  of  obtaining  thrust  and  aerodynamic  drag  from  reduction  of 
camera  data  with  the  aid  of  high  speed  computers  are  given.  Determi¬ 
nation  of  various  performance  parameters  needed  for  the  construction 
of  firing  tables  is  explained;  some  typical  results  are  presented. 
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TABLE  OF  SYMBOLS 


Term 

Definition 

Unit 

A 

e 

Area  of  jet  exit 

ft3 

AZ 

Azimuth  of  line  of  fire  (clockwise  from  North) 

deg 

B  (t) 

Transverse  moment  of  inertia  at  time  t 

lbm-ft3 

C 

Ballistic  coefficient 

lbm/ina 

d 

Reference  diameter  of  missile 

ft 

E 
— * 

Position  of  the  missile  with  respect  to 
spherical  Earth  surface 

ft 

F 

Deceleration  due  to  friction 

ft/sec  3 

£ 

Acceleration  vector  due  to  gravity 

ft/sec3 

go 

Acceleration  due  to  gravity  (surface) 

ft/sec3 

gc 

Constant  used  in  the  conversion  of 
lbf-sec3  /ft  to  lbm 

lbm-ft/  lbf-sec3 

h 

-♦ 

Angular  momentum  divided  by  B  (t) 

rad/ sec 

h*i 

ha  at  end  of  launch  (t^  ) 

rad/ sec 

hai 

h3  at  end  of  launch  (t^ ) 

rad/sec 

H 

Total  angular  momentum 

lbm-ft3  -rad/sec 

• 

H 

-4 

Rate  of  change  of  H 

lbm-ft3  -rad/sec 

*SP 

Specific  impulse  per  unit  of  fuel  mass 

lbf-sec/  lbm 

!st 

Total  impulse  at  standard  conditions 

Effective  drag  force  coefficient 

lbf-sec 

kd 

V 

Drag  force  coefficient  during  thrust-on 

db 

\r 

Zero  yaw  drag  force  coefficient 

J\_ 

D0 

kd 

a 

kh 

Yaw  drag  force  coefficient 

Damping  moment  coefficient 

l/rad3 
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TABLE  OF  SYMBOLS  (Contd) 

Term  Definition  Unit 

Lift  force  coefficient  ■  ■— — — 

Overturning  moment  coefficient  —  ■ 

K  Pitching  force  coefficient  » . 


L  Latitude  of  launch  point  deg 

m  Mass  of  missile  lbm 

m  (t)  Mass  of  missile  at  time  t  lbm 

m^  Mass  of  missile  at  end  of  thrust-on 

BST  phase  at  standard  conditions  lbm 

m^  Mass  of  fuel  lbm 

m  Mass  of  fuel  at  standard  conditions  lbm 

fST 

m  (t)  Rate  of  change  of  mass  at  time  t  lbm/sec 

P^  Static  atmospheric  pressure  lbf/ft3 

P^  Jet  pressure  at  nozzle  exit  lbf/ft 3 

PE  Probable  error  — 

PT  Propellant  temperature  #F 

QE  Quadrant  elevation  r/ 

r  Distance  between  center  of  Earth  and 

projectile  ft 

r  Distance  from  CG  of  missile  to  nozzle  exit  ft 

e 

r  Distance  from  CG  of  missile  to  nozzle 

*  throat  ft 

R  Radius  of  Earth  ft 

t  Time  sec 

tg  Motor  burnout  time  sec 

t__ ,  Time  from  motor  ignition  until  first  motion  sec 

BM 

*BMST  Time  from  motor  ignition  until  first  motion 

at  standard  conditions  sec 
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TABLE  OF  SYMBOLS  (Contd) 


Term 


tBST 

T 

T* 


tr“> 

u 

-* 

u 

—4 

V 

-4 

w 

-♦ 

X 

-4 

X 

-4 

X 

-4 


a 

A 

-4 

P 

0 

n 


X 


Definition 

Motor  burnout  time  at  standard  conditions 
Thrust 

Effective  thrust 
Thrust  factor 

Thrust  produced  by  motor  at  time  t 

Velocity  of  missile  with  respect  to  ground 
Acceleration  of  missile  with  respect  to  ground 
Velocity  of  missile  with  respect  to  air 
Velocity  of  air  with  respect  to  ground 
Unit  vector  along  longitudinal  axis  of  missile 
Rate  of  change  of  x 

Position  of  missile  with  respect  to  ground- 
fixed  coordinate  system 

Yaw  of  missile 

Coriolis  acceleration  due  to  rotation  of  Earth 

Air  density  (varies  with  altitude) 

Angle  of  elevation  of  launch 

Angular  velocity  of  the  Earth 

Dot  -  when  superscripted,  eg.  u,  denotes  first 
derivative  as  a  funtion  of  time  -  when  appearing 
between  vectors,  eg  v  .  x  ,  denotes  a  dot 
product 

Denotes  cross  product  of  vectors 


Unit 

sec 

lbf 

lbf 


lbf 

ft/sec 

ft/sec3 

ft/sec 

ft/sec 


rad/ sec 
ft 

deg 

ft/  sec3 

lbm  /ft3 
deg 

rad/sec 
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PAGES  NOT  FI1M3D  ARE  BLANK 
INTRODUCTION 

The  purpose  of  this  report  is  to  describe  the  methods  used  by  the 
Ballistic  Research  Laboratories  in  determining  performance  parameters 
for  fin-stabilized,  free-flight  missiles  for  use  in  generating  firing  tables. 

The  work  presented  herein  is  derived  from  a  simplified  version  of  a 
general  rigid  body  mathematical  model  described  in  BRL  Report  No.  1244.  1  * 

FREE-FLIGHT  MISSILE  TEST  PROGRAM 

A  firing  program  for  production  missiles  is  designed  to  determine 
the  performance  of  a  free-flight  missile  system.  This  program,  or 
matrix,  outlines  the  specific  conditions  of  quadrant  elevation  (QE)  and 
propellant  temperature  (PT)  under  which  the  rounds  are  fired.  A 
typical  matrix  defining  a  test  program  is  shown  below. 


Number  of  Missiles 

PT 

QE  (r*0  -40  20 

(”F) 

77 

120 

100 

4 

4 

4 

4 

200 

4 

4 

4 

4 

400 

4 

4 

4 

4 

800 

4 

4 

4  ” 

4 

This  array  permits  determination  of  motor  performance  as  a  function  of 
propellant  temperature,  and  launcher-missile  performance  as  a  function 
of  quadrant  elevation.  The  missiles  are  fired  singly  from  at  least  two 
launchers. 

The  following  data  are  collected  for  each  round. 

1.  Missile  pre-flight  measurements 

a.  Weight  of  components 

b.  Center  of  gravity,  with  and  without  propellant 

c.  Moment  of  inertia,  with  and  without  propellant 

Superscript  numbers  denote  references  which  may  be  found  on  page  29. 
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d.  Propellant  temperature 

e.  Fuze  setting 

f.  Warhead  type 

2.  Missile  in-flight  measurements 

a.  Positions,  velocities,  and  accelerations  from  camera  data 

(1)  From  ignition  to  clear  of  launcher  at  100  frames  per  second 

(2)  From  launcher  to  point  beyond  maximum  velocity  at  30 
frames  per  second 

(3)  From  maximum  velocity  to  terminal  event  at  5  frames  per 
second 

b.  Roll  data 

3.  Missile  post-flight  measurements 

a.  Position  of  impact  (range,  height  and  deflection) 

b.  Time  of  impact 

4.  Launcher  measurements 

a.  Tactical  emplacement,  before  and  after  firing 

b.  Quadrant  elevation,  before  and  after  firing 

c.  Azimuth  of  fire,  before  and  after  firing 

d.  Survey  of  position,  before  and  after  firing 

e.  Type  and  serial  number 

5.  Meteorological  measurements 

a.  Winds 

(1)  Recorded  at  five-second  intervals  at  the  launch  site  before, 
during,  and  after  firing  from  0  to  100  feet  above  the  launcher 

(2)  By  balloon  at  a  location  and  time  as  close  to  the  launch  site 
and  firing  time  as  possible 

(a)  From  0  to  1000  feet  above  the  launcher  recorded  at 
thirty-second  intervals  during  ascent 

(b)  From  surface  to  a  height  greater  than  the  expected 
maximum  ordinate  of  the  missile  recorded  at  sixty-second 
intervals  during  ascent 

b.  Atmospheric  data  (pressure,  temperature  and  density) 

(1)  At  the  launch  site  at  the  time  of  fire 
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(2)  By  balloon  at  a  location  and  time  as  close  to  the  launch  site 
and  firing  time  as  possible  from  surface  to  a  height  greater 
than  the  expected  maximum  ordinate  of  the  missile  recorded 
at  sixty-second  intervals  during  ascent. 

REDUCTION  OF  FLIGHT  TEST  DATA  FOR  DRAG  COEFFICIENT  (KD) 

AND  THRUST (T)  FORCES 

The  drag  force  and  ultimately  the  drag  coefficient  and  the  thrust 
forces  are  determined  by  an  inverse  solution  of  the  point-mass 
differential  equations  of  motion.  The  drag  force  is  obtained  from  the 
free-flight  portion  of  the  trajectory,  and  the  thrust  force  is  found 
from  the  boost  phase.  Two  assumptions  are  made  to  carry  out  the 
latter  phase. 

1.  The  drag  force  during  boost  is  small  in  comparison  to  the  thrust 
force  and  is  the  same  as  that  found  for  free-flight  with  estimated  changes 
to  the  drag  force  to  account  for  the  presence  of  the  jet. 

2.  An  initial  thrust  curve  is  assumed  along  with  the  following  equation 
to  determine  the  mass  flow  rate. 


(1.1)  m  (t) 
where  m  (t) 

trU) 


-tr  (t)/igP, 

Rate  of  change  of  mass  at  time  t  , 
Thrust  produced  by  motor  at  time  t  and 

Specific  impulse  per  unit  of  fuel  mass . 


The  equations  used  for  the  determination  of  the  drag  and  thrust 
forces  are 

(1.2)  Kd  =  -  C  (ux  -  wx)  (uj,  +  gx  +  Ai )  +  u3  (ua  +  g3  +  Aa) 
+  (u3  -  w3  )  (U3  +  g3  +  A3 )  ]  m/ p  da  va  and 
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(1.3)  T  =  {  C  (ui  -  wx)  (ux  +  gi  +  Ai)  +  ua  (ua  +  ga  +  A3)  +  (u3  -  w3) 


(u3  +  g3  +  As)  ]  m  (t)/g  c  v  +  (p  da  va/gc  )  kd 


The  approximation  of  the  force  of  gravity  is 
(1.4)  jj|  s  ( “go  Ra/r3) 

where  r  =  [X?  +  (  Xa  +  R  )a  +  X3a  ]* , 

go  =  Value  of  gravity  at  point  of  launch, 
r  =  Distance  between  center  of  Earth  and  projectile, 
R  =  Radius  of  Earth,  and 

X  =  Position  of  missile  with  respect  to  ground-fixed 
coordinate  system  . 

The  Coriolis  acceleration  due  to  rotation  of  the  Earth  is 


Xa  +  R 
X3 


(1.5) 


A  i  ua  —  Aa  u3 
A  i  u^  4  A3  u3  • 
A  a  nl  —  A3  ua 


For  the  Northern  Hemisphere  the  A's  are  defined  by  the  following 
equations  [  for  the  Southern  Hemisphere  replace  L  by  (-L)  ]  : 


Aj_  =  ZOcos  L  sin  AZ, 
A a  =  2  Osin  L  and 


A3  =  2  ftcos  L  cos  AZ , 


where  AZ  =  Azimuth  of  line  of  fire  (clockwise  from  North)  , 
L.  =  Latitude  of  launch  point  and 

Q  =  Angular  velocity  of  the  Earth. 
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The  mass  of  the  individual  missile,  the  atmospheric  conditions 
at  time  of  firing  and  the  reduced  camera  data  provide  the  necessary 
inputs  to  solve  Equations  (1.2)  and  (1.3)  for  drag  and  thrust.  Examples 
of  the  results  of  this  method  of  reduction  are  presented  as  Figures  1  and  2. 

The  thrust-time  history  of  a  rocket  motor  as  a  function  of  propellant 
temperature  is  obtained  by  assuming  a  linear  transformation  of  the 
standard  propellant  temperature  thrust  curve  in  time. 


Let 

(1.6) 


where 


T 


/  tBST  *BMSt\  + 
'  tB  _tBM  / 


f  [  ]  indicates  "  a  function  of  " 


tBMST 
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The  thrust  of  the  rocket  motor  as  a  function  of  time  may  be  linearly 
transformed,  maintaining  the  same  total  impulse,  by  the  following  equation: 


(1.7) 


TR  =  mf 


/  I_SP  \  /  tBST  ~*BMST  \ 

\  !ST  /  \  tB  -tBM  / 


T. 


The  force  due  to  jet  pressure  at  the  rocket  nozzle  is  given  by  P  A  . 

C  6 


The  static  atmospheric  pressure  produces  a  force  equal  to  -P  A  . 

cl  6 

The  total  thrust  contributed  by  pressure  is  then  given  by  (P  -P  )  A  . 

0  EL  G 

The  effective  thrust  of  the  rocket  motor  is  the  algebraic  sum  of  the 
components  considered,  and  is  given  by 


(1.8) 


T*  -  t  m  fISP)ftBST  BMST  \ 


+  (P  -P)  A 
e  a  e 


FLIGHT  SIMULATION  OF  FREE-FLIGHT  MISSILES 


The  thrust  curve  as  a  function  of  time  and  the  drag  force  coefficient 
curve  as  a  function  of  Mach  number  are  combined  with  other  aerodynamic 
characteristics  (obtained  from  wind  tunnel  testing,  inflight  spark  range 
measurements  and  estimates  based  on  flow  calculations)  for  use  in  the 
following  phases  which  are  explained  in  detail  in  BRL  Report  No.  1244:1 

1.  Launching  phase 


(2.1)  u  = 


^  c  T  cos  0  -pda  K  vs  cos  0 
m  (t)  m  (t)  B 

-  go  sin  0  cos  0  -F  cos3  0 


8  c  sin  0  -pd3  K  v3  sin  0 
m  (t)  m  (t)  B 

-  go  sin3  0  -F  cos  p  sin  0 
0 
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2.  Boost  phase  (rigid  projectile) 

m  (t) 


(2-2)  2  -  -  <KDll  +  kd  °',)vX  +  Sl7)  klCvS  5-<S* 

B  OL 


+-£4 


m(t)  B(t) 


Ksv(xXH)  +  *  +  4  + 


g  T* 
m  (t) 


(2.3)  U  =  pd3  Km  v(vX*j  +  Kh  v  [(H.  x)  x  -  H] 


(m  (t)  r  r  ) 

+  - p  -  [H  -  (H  .  x)  x  ]  and 

B  (t)  -*  -♦ 


(2.4)  x  =  (HXx)/B(t). 

3.  Free-flight  phase  (point  mass) 

(2.5)  u  =  -  (pm  /144  m  C)  v  v  +  £  + 

BST 


(2. 6)  X  =  u  dt  and 

X! 

Xa  +  R  -  (R3  -  X?^  . 

X3 

In  the  interest  of  economy,  flight  simulations  are  carried  out  by  two 
distinct,  but  combined,  phases.  The  boost  portion  is  simulated  by  the  use 
of  the  rigid-body  model  mentioned  above,  and  the  point-mass  model  is 
used  in  the  free-flight  phase.  Since  any  initial  yawing  motion  of  a  fin-stabi¬ 
lized,  free-flight  missile  is  essentially  damped  out  by  motor  burnout  time, 
no  degradation  of  result  occurs  from  j simulation  of  the  free-flight  phase 
with  a  point-mass  model. 


(2.7)  E  = 
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Individual  missile  weights  and  measures  and  the  atmospheric  condi¬ 
tions  at  the  time  of  firing  are  used  by  the  simulator  with  certain  chosen 
variables  to  match  the  flight  of  the  missile  at  two  key  points  along  the 
trajectory.  These  key  points  are  obtained  by  the  camera  instrumentation 
and  are 

1.  Velocity  vector  at  a  time  shortly  after  motor  burnout  and 

2.  Range  and  height  of  terminal  event. 

The  velocity  vector  is  matched  iteratively  by 

1.  A  multiplier  (T  )  on  thrust  applied  to  the  thrust  curve  to  obtain 
the  magnitude  of  the  velocity  vector  and 

2.  Initial  angular  velocities  in  the  vertical  (h3^)  and  horizontal  (ha^) 

planes  at  the  time  of  launch  to  match  the  total  angle  of  the  velocity 
vector.  These  are  merely  pseudo  rates  to  enable  the  simulator  to 
match  the  missile  attitude. 

The  range  is  matched  iteratively  by  the  simulator  through  the  use  of  an 
inverse  multiplier  (C)  on  the  drag  curve. 

Some  additional  comparisons  are  made  to  determine  exactness  of  the 
simulation.  These  are 

1.  Velocity  at  end  of  launch, 

2.  Position  of  missile  at  time  of  match  of  velocity  vector  , 

3.  Position  of  missile  and  velocity  vector  at  a  point  approximately 
2  seconds  beyond  time  of  match  of  velocity  vector  » 

4.  Position  of  missile  and  velocity  vector  near  maximum  ordinate,  and 

5.  Deflection  and  time  of  flight  at  terminal  event . 

When  these  comparisons  are  made  to  the  degree  of  accuracy 
desired,  the  following  performance  parameters  for  a  free-flight  missile 
are  obtained: 

1.  Effective  drag  force  coefficient  as  a  function  of  Mach  number  as 
determined  from  the  drag  reduction.  (Fig.  3) 

2.  Thrust  as  a  function  of  time  as  determined  from  the  thrust  re¬ 
duction.  (Fig.  4) 
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3.  Burning  time  of  motor  as  a  function  of  propellant  temperature 
as  determined  from  the  thrust  reduction.  (Fig.  5) 

4.  Thrust  factor  (T  )  as  a  function  of  propellant  temperature. 

(Fig.  6)  This  represents  the  effect  of  propellant  temperature  on 
the  total  impulse  of  the  rocket  motor;  therefore,  the  observed 
probable  error  in  thrust  factor  is  representative  of  the  probable 
error  in  impulse. 

5.  Initial  angular  velocities  (h and  h3^)  as  functions  of  quadrant 
elevation.  (Figs.  7  and  8) 

6.  Ballistic  coefficient  (C)  as  a  function  of  quadrant  elevation.  (Fig.  9) 

On  Figures  6  through  9  the  solid  line  indicates  the  mean  performance 
values  of  the  particular  parameter;  whereas, the  dashed  lines  present  an 
idea  of  the  dispersion  of  that  parameter. 

CONCLUSION 

The  methods  described  above  have  been  used  successfully  to  describe 
the  Littlejohn  and  Honest  John  flight  performance  and  in  preparing  adequate 
firing  tables  for  those  weapon  systems. 


ROBERT  F.  LIESKE 


JOSEPH  W.  KOCHENDERFER 


Motor  Burnout  Time  (tB)vs  Propellant  Temperature  (PT) 


Figure 


Angular  Velocity  at  End  of  Launch- Vertical  (ha^  )  vs  Quadrant  Elevation  (QE) 


s 
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